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Abstract. In view of the extensive evidence of a tight inter-relationship between 
spheroidal galaxies (and galactic bulges) and massive black holes hosted at their cen¬ 
ters, a consistent model must deal jointly with the evolution of the two components. 
We describe one viable model, which successfully accounts for the local luminosity 
function of spheroidal galaxies, their photometric and chemical properties, deep galaxy 
counts in different wavebands, including those in the (sub)-mm region which proved to 
be critical for current semi-analytic models stemming from the standard hierarchical 
clustering picture, clustering properties of SCUBA galaxies, of EROs, and of LBGs, 
as well as for the local mass function of massive black holes and for quasar evolution. 
Predictions that can be tested by surveys carried out by SIRTF are presented. 


1 Introduction 


The hierarchical clustering model with a scale invariant spectrum of density 
perturbations in a Cold Dark Matter (CDM) dominated universe has proven to 
be remarkably successful in matching the observed large-scale structure as well 
as a broad variety of properties of galaxies of different morphological types (e.g. 
[|6|Jl^). However, serious shortcomings of this scenario have also become evident 
in recent years. The critical point can be traced back to the relatively large 
amount of power on small scales predicted by this model which would imply 
far more dwarf galaxies or substructure clumps within galactic and cluster mass 
halos than are observed (the so-called “small-scale crisis” ([pOjM,^,^), unless 
star formation in small objects is strongly suppressed (or the small scale power 
is reduced by modifying the standard model). 

At the other extreme of the galaxy mass function we have another strong 
discrepancy with model predictions, that we might call “the massive galaxy 
crisis”: even the best semi-analytic models (||^,^) hinging upon the standard 
picture for structure formation in the framework of the hierarchical clustering 
paradigm, fall short by a substantial factor (up to about 10) to account for the 
(sub)-mm (SCUBA and MAMBO) counts of galaxies, most of which are prob¬ 
ably massive objects undergoing a very intense star-burst (with star formation 
rates ^ lOOOMQyr"^) at z > 2 (see, e.g. [|^). Recent optical data confirm that 
most massive ellipticals were already in place and (almost) passively evolving 
up to z ~ 1-1.5, implying that they were fully assembled by z ^ 2.5, although 
the issue is still somewhat controversial These data are 
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more consistent with the traditional “monolithic” approach whereby giant el¬ 
lipticals formed most of their stars in a single gigantic starburst at substantial 
redshifts, and underwent essentially passive evolution thereafter. 

On the contrary, in the canonical hierarchical clustering paradigm the small¬ 
est objects collapse first and most star formation occurs, at relatively low rates, 
within relatively small proto-galaxies, that later merged to form larger galaxies. 
Thus, the expected number of galaxies with very intense star formation is far 
less than detected in SCUBA and MAMBO surveys and the surface density of 
massive evolved ellipticals at 2 ^ 1 is also smaller than observed. The “mono¬ 
lithic” approach, however, is inadequate to the extent that it cannot be fitted in a 
consistent scenario for structure formation from primordial density fluctuations. 


2 Relationships between quasar and galaxy evolution 


The above difficulties, affecting even the best current recipes, may indicate that 
new ingredients need to be taken into account. A key new ingredient may be 
the mutual feedback between formation and evolution of spheroidal galaxies and 
of active nuclei residing at their centers ( |p^ , p^ , ^l|J^Jl^ , |l^,|l6y3lj^ , ^ ). In this 
framework, |l^ elaborated the following scheme (see also lliy^,|3i^): 


Feed-back effects, from supernova explosions and from active nuclei (note 
that supernova feedback alone falls short of solving the dearth of dwarf 
galaxies, [^, but photo-ionization by the UV background re-ionizing the 
inter-galactic medium (IGM) could do the job |M ) delay the collapse of 
baryons in smaller clumps while large ellipticals form their stars as soon as 
their potential wells are in place; the canonical hierarchical CDM scheme - 
small clumps collapse first - is therefore reversed for baryons. 

Large spheroidal galaxies therefore undergo a phase of high (sub)-mm lumi¬ 
nosity. 

At the same time, the central black-hole (BH) grows by accretion and the 
quasar luminosity increases; when it reaches a high enough value, its action 
(ionization and heating of the gas) stops the star formation and eventually 
expels the residual gas. This explains the observed correlation between BH 
and host spheroidal masses (see [ ^3|jl^ ). The same mechanism distributes 
in the IGM a substantial fraction of metals and may pre-heat the IGM. The 
onset of quasar activity (and the corresponding squelching of star forma¬ 
tion) occurs earlier for more massive objects. The duration of the star-burst 
increases with decreasing mass from ~ 0.5 to 2Gyr. 

This implies that the star-formation activity of the most massive galaxies 
quickly declines for z 3, i.e. the redshift distribution of SGUBA/MAMBO 
galaxies should peak at z 3, just before quasars reach their maximum 


luminosity (at z ~ 2.5). This: 

— explains why very luminous quasars are more easily detected at (sub)- 
mm wavelengths for z ^ 2.5 [^. The latter authors argue that a large 
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fraction of the observed (sub)-mm emission is powered by a starburst, 
as expected from this model 

— implies an extremely steep (essentially exponential) decline of (sub)-mm 
counts at bright fluxes, as indicated by recent data. 

• A “quasar phase” follows, lasting 10^-10® yrs. 

• A long phase of passive evolution of galaxies ensues, with their colors be¬ 
coming rapidly very red [Extremely Red Object (ERO) phase]. 

• Intermediate- and low-mass spheroids have lower Star Formation Rates and 
less extreme optical depths. They show up as Lyman-Break Galaxies (LBGs). 

• Therefore, in this scenario, large ellipticals evolve essentially as in the ’’mono¬ 
lithic” scenario, yet in the framework of the standard hierarchical clustering 
picture. 

The various aspects and implications of this compound scheme have been ad¬ 
dressed by our group in a series of papers. |5ll estimated the mass function of 
quiescent BHs at the centers of local galaxies. They found a dichotomy between 
large and small BHs: the former are hosted in elliptical galaxies and tend to be 
not obscured, while the latter, found in the bulges of spiral galaxies, can be reac¬ 
tivated by interactions and are frequently obscured. Their results are consistent 
with the high luminosity quasar activity occurring as a single short-lived event at 
relatively high redshift with luminosity close to the Eddington limit. Later, lower 
luminosity, nuclear activity may be due to re-activation of the central BH, due, 
e.g., to interactions. analyzed the evolution of the quasar luminosity func¬ 
tion in the framework of a model in which the spheroidal galaxies (and bulges 
of later-type galaxies) and the BHs at their centers form and evolve in parallel. 
Adopting the standard |Q formation rate of dark matter halos, they found that 
consistency with the data is achieved if, and only if, BHs in bigger galactic halos 
form earlier. In other words, the interval between the onset of star formation 
and the peak of quasar luminosity is shorter for the more massive objects. Some 
of the other, most recent, results are briefly described in the next section. 

3 Some recent results and predictions 

3.1 Counts at (sub)-mm wavelengths 

The (sub)-mm counts are expected to be very steep because of the combined ef¬ 
fect of the strong cosmological evolution of dust emission in spheroidal galaxies 
and of the strongly negative K-correction (the dust emission spectrum steeply 
rises with increasing frequency). The model by has extreme properties in 
this respect: above several mJy its 850 /im counts reflect the high-mass exponen¬ 
tial decline of the mass function of dark halos. In this model, SGUBA/MAMBO 
galaxies correspond to the phase when massive spheroids formed most of their 
stars at z ^ 2.5; such objects essentially disappear at lower redshifts. On the 
contrary, the counts predicted by alternative models (which are essentially phe¬ 
nomenological [p| J^J^ , ^ , p^ ), while steep, still have a power law shape, and the 
redshift distribution has an extensive low-z tail. 
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Fig. 1. Integral source counts at 850/rm (left panel) and at 1.4 mm (right panel) pre¬ 
dicted by the model by [^] compared with observations. The dotted, dashed and dot- 
dashed lines show the contributions of starburst, spiral, and forming elliptical galaxies, 
respectively. The long-dashed line (shown only at 1.4 mm) gives the counts of radio 
sources, after |)7I|. T he solid line shows the total counts. References for the data points 
can be found inT^] 


As illustrated by Fig. 1, the recent relatively large area surveys are 

indeed suggestive of an exponential decline of the 850 /rm counts above several 
mjy. Further evidence in this direction comes from MAMBO surveys at 1.2 mm 
§; see ||). 


3.2 Leasing effects on (sub)-mni counts 

A direct consequence of the extreme steepness of the (sub)-mm counts predicted 
by this model is that their bright tail is strongly affected by gravitational leasing 

In fact, although the probability of strong leasing is very small, it has a 
power-law tail (p(A) oc A“^) extending up to large values of the magnification A. 
Thus, if counts are steep enough, the fraction of leased sources at bright fluxes 
may be large. 

and iQ find that, according to the model by |^, essentially all proto- 
spheroidal galaxies brighter than Sgso^m — 60-70 mJy are gravitationally leased. 
Allowing for the other populations of sources contributing to the bright mm/sub- 
mm counts, they find that the fraction of gravitationally leased sources may be 
~ 40% at fluxes slightly below S'sso/im = 100 mJy. If so, large area surveys such as 
those to be carried out by Planck/HFI or by forthcoming balloon experiments 
like BLAST and ELISA will probe the large scale distribution of the peaks of 
the primordial density field. For comparison, the maximum fraction of leased 
sources predicted by current phenomenological models is 5%. 
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3.3 Clustering of SCUBA/MAMBO galaxies 

Since, in the above scheme, SCUBA galaxies correspond to the rare, massive 
density peaks at high-z, they are expected to be highly biased tracers of the dark 
matter distribution, and therefore to be strongly clustered. The recent analyses 
by and [Q have shown that the implied angular correlation function is 
indeed consistent with the results by 0 on the clustering of LBGs, and by 
[HI and @ on clustering of SCUBA galaxies. The estimated w{9) for bright 
SCUBA galaxies HI and for EROs ^ indicate dark halo masses 10^^ Mq for 
these sources; the w{6) for LBGs Hi indicate typical masses at least 10 times 
lower (see also 0). 


3.4 Chemical and photometric evolution 

The problem has been addressed by [0. The short duration of the star-burst 
in spheroidal galaxies, and its occurrence at z 3, assumed by the model, are 
consistent with the supra-solar Mg/Fe ratio in ellipticals, with the tightness 
of their fundamental plane and of their color-tr relation. But the model also 
compares successfully with a broad variety of detailed observational data: the 
correlation between the abundances of Fe and those of “enhanced” elements 
(C, N, O, Ne, Mg, Si, S), the correlation between Fe and Mg 2 , the template 
Spectral Energy Distribution of local ellipticals, the color (V — A)-magnitude 
(My) and the {J — K)-{V — K) relations. The key ingredients allowing to match 
the observational data, and in particular the relationships among abundances of 
the various elements are, on one side, the decrease of the duration of the star- 
burst with increasing galactic mass, and, on the other side, the fact that the 
effect of stellar feedback is higher in less massive objects, so that the amount 
of material enriched by SNII that can be retained is lower. The main difference 
with the majority of previous chemical evolution models is that star formation 
is stopped not by supernova-driven winds (which also fall short of providing 
enough preheating of the proto-cluster medium, 00,01), but by the energy 
injected by the active nucleus. 

Note that this model also accounts for the roughly solar, or supra-solar, 
metallicity of the quasar broad-line regions and of supra-solar N/C and Feja 
ratios in the more luminous objects observed out to z > 4 |0| and even up to 
z ^ 6 [0. This implies that the medium must have been chemically enriched 
before the quasars became observable and on short timescales (< 1 Gyr, at least 
at the highest redshifts). 

Another interesting implication of the model, also discussed by [0, is the 
dependence of the ratio Mdarkmatter/Mgtars on Mgtars- This ratio has a minimum 
value ~ 20 for Mgtars — 10^°-10^^ Mq, in good agreement with the results by 
[0, and increases at smaller masses (in keeping with the findings by [0) due 
to the increasing fraction of gas expelled by supernova feedback, as well as at 
higher masses, due to the increase of the cooling time of the gas in the outer 
regions of the galactic halo. 
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Fig. 2. Integral source counts (left-hand panel) and redshift distribution (right-hand 

S iel) for a flux limit of 30/rJy (right panel) at 24 fim predicted by the model by 
]. In the left-hand panel, the dotted, dashed and dot-dashed lines show the con¬ 
tributions of starburst, spiral, and forming elliptical galaxies, respectively, while the 
three-dots/dashed line shows the total counts of ellipticals, including also those where 
the star-formation has ended. In the right-hand panel, the solid and the dotted lines 
show the redshift distributions of ellipticals during and after the star-formation phase, 
respectively 


3.5 Predictions for SIRTF surveys 


SIRTF surveys have the potential of providing further tests of the model. In par¬ 
ticular, the GOODS (http://www.stsci.edu/science/goodg) 24survey should 
reach the confusion limit at 30--100 ^Jy. According to the model, about 50% of 
the detected galaxies should be spheroidal galaxies forming their stars a,t z ^ 2. 


About 400-600 such objects are expected over an area of 0.1 square degree (see 
Fig. 2). Their redshift distribution is predicted to peak at z slightly above 2, 
with a significant tail extending up to z 3. 
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